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For two-dimensional lattices in a tight-binding description, the intrinsic spin-orbit coupling, acting as a com- 
plex next-nearest-neighbor hopping, opens gaps that exhibit the quantum spin Hall effect. In this paper, we 
study the effect of a real next-nearest-neighbor hopping term on the band structure of several Dirac systems. In 
our model, the spin is conserved, which allows us to analyze the spin Chern numbers. We show that in the Lieb 
and kagome lattices, variation of the amplitude of the real next-nearest-neighbor hopping term drives interesting 
topological phase transitions. These transitions may be experimentally realized in optical lattices under shaking, 
when the ratio between the nearest- and next-nearest-neighbor hopping parameters can be tuned to any possible 
value. Finally, we show that in the honeycomb lattice, next-nearest-neighbor hopping only drives topological 
phase transitions in the presence of a magnetic field, leading to the conjecture that these transitions can only 
occur in multi-gap systems. 
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I. INTRODUCTION 

The study of topological states of matter has received a 
large increase of interest since the theoretical formulation 1 2 
and experimental observation 3 of the quantum spin Hall ef- 
fect. Besides the quantum Hall and quantum spin Hall 
effects in two dimensions, several other topological states 
of matter have been studied, e.g., topological insulators in 
three dimensions 45 and topological superconductors^ and 
they have been classified in a systematic fashion. 7 This clas- 
sification is based on the notion of topological invariants; in 
the quantum Hall and spin Hall systems, these invariants are 
closely related to the quantized Hall and spin Hall conductiv- 
ity, respectively. These conductivities can be observed when 
the bulk is in an insulating phase; they are edge properties that 
exist inside the bulk gaps. Remarkably, these quantities are 
closely related to the topological invariants of the bulk bands 
known as the Chern numbers, a property known as the bulk- 
boundary correspondence. Importantly, these Chern numbers 
make it possible to classify the topological phases of the sys- 
tem, without having to study the physics on the edges. 

Phase transitions between these topological phases are as 
appealing as the phases themselves. Topological phase tran- 
sitions are characterized by a change of the topological in- 
variants that classify the system. Because of the topological 
protection, they can occur only when a gap closes. In that 
case, the Chern numbers of the energy bands change. 4 Several 
mechanisms can give rise to these transitions. For instance, 
in the model by Kane and Mele, 1 it has been shown that by 
increasing the Rashba spin-orbit coupling one can close the 
quantum spin Hall gap and open a trivial one. Other types of 
transitions, between different chiral states, occur in the Hofs- 
tadter model, and are driven by variation of the magnetic flux 
density. 8 More recently, it has been shown that the intrinsic 
spin-orbit coupling also drives phase transitions in the honey- 
comb lattice in the presence of magnetic flux I2U21 In this work, 
we show that in several two-dimensional lattices, topologi- 
cal phase transitions can be driven by tuning the strength of 
a real next-nearest-neighbor (NNN) hopping term. Although 
in tight-binding models this term resembles the intrinsic spin- 



orbit (ISO) coupling (which is a purely imaginary NNN hop- 
ping term), it is often neglected, despite the fact that in some 
systems it may be much stronger than the ISO coupling. 

Experiments incorporating ultracold fermionic atoms in op- 
tical lattices are ideal candidates for the realization of topolog- 
ical phase transitions, since the strengths of the different hop- 
ping parameters can be tuned relatively easily. In the mostly 
realized square optical lattices, the potential in the x and y 
directions is usually separable and the NNN hopping (along 
the diagonal) is zero. For non-separable optical lattices, the 
magnitude of the NNN hopping t' is small compared to that 
of the nearest-neighbor (NN) hopping t. However, it has been 
recently shown that by shaking the optical lattice the effec- 
tive ratio t' cS /teg can be varied in the entire range from to 
oo. This occurs because, in a certain regime of parameters, 
the shaking leads to a renormalization of the hopping, which 
becomes multiplied by a Bessel function, i.e., t e g = tJo{K), 
where K — Vq/Hu! is the shaking parameter, given by the ra- 
tio between the amplitude Vo and the frequency u) of the peri- 
odic shaking. 11 Moreover, the NNN may get renormalized by 
a Bessel function with a different argument, thus leading to a 
full range of possible values for the effective ratio t' eS /t e ffP^ 

Here, we focus on tight-binding models, where we assume 
that all particles loaded onto the lattice are in their respec- 
tive ground states, and we study topological phase transitions 
driven by tuning the mtiS^t' jt in the presence of ISO cou- 
pling. In this aspect, we study several lattice geometries. The 
honeycomb lattice (e.g., of graphene) is one of the most fa- 
mous systems exhibiting a Dirac dispersion,^ where ISO cou- 
pling leads to the quantum spin Hall phase. However, the fact 
that there are only two (spin-degenerate) bands, and hence 
only one gap, greatly reduces the possible types of topolog- 
ical phase transitions that could occur. In the kagome and 
Lieb lattices, the spectrum consists of three bands, and hence 
two gaps, increasing the number of possible topological phase 
transitions. In this work, we discuss these three lattices and 
compare the topological phases and topological phase transi- 
tions that they exhibit. 

The general tight-binding Hamiltonian, that applies to all 
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FIG. 1. (Color online) (a) The geometry of the Lieb lattice. The 
square unit cell with lattice constant a = 1 is indicated, together with 
the three sublattices (A, B, C). An example of a NN (NNN) hopping 
process t (t') is indicated by the red arrows, (b) First Brillouin zone 
of the Lieb lattice. The high symmetry points F, X, Y and M are 
given. The red dashed lines indicate the path TXMT along which 
the dispersions are displayed in the subsequent figures. 

lattices presented in this work, reads 

H = i?NN + ^NNN ~T" Hiso, (1) 

with 

^NNN = ~t' 4,<A>) (3) 

HlSO = +«Ai SO ^ Si,A e ij • f<T(J')Sj,<r'- (4) 

Here, Si_ a (s\ a ) is the annihilation (creation) operator for a 
particle on site i and with spin a. The first two terms i?NN 
and -£/nnn describe NN and NNN hopping, respectively. The 
third term i/iso is the ISO coupling; in this term, we have 
defined the unit vector = (d.jj. x djy)/|djfc x | in terms 
of the bond vectors d^ and dy, that connects the sites i and 
j via the unique intermediate site k. Finally, cr = [a x , cr y ,a z ) 
is the vector of Pauli matrices. 

This paper is organized as follows. In Section [II] we apply 
the tight-binding model to the Lieb lattice, study the spectrum 
and analyze the topological phase transitions using the Chern 
numbers of the bands. In Section [Ill] w e repeat this discussion 
for the kagome lattice. In Section jTV] we study a honeycomb 
lattice and show that contrarily to the other two examples, no 
topological phase transitions driven by the real NNN hopping 
can occur, unless a magnetic field is applied perpendicularly 
to the lattice. Finally, we discuss the possibilities for experi- 
mental realization of this model and we conclude in SectionM 

H. THE LIEB LATTICE 

A. Intrinsic spin-orbit coupling and next-nearest-neighbor 
hopping 

First, we investigate various topological phase transitions 
on the Lieb lattice. In Fig.[TJa), we show the lattice geometry, 
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FIG. 2. (Color online) (a) The upper and lower bands touch the mid- 
dle flat band at the M -point if Aiso = t' — 0. (b) A gap appears if 
the ISO coupling is turned on (Aiso 7^ and t' = 0). 



with three sites A, B, and C per unit cell, and lattice con- 
stant a = 1. The high-symmetry points of the first Brillouin 
zone (1BZ) are shown in Fig. [TJb). Both in real space and 
in reciprocal space, the system satisfies C4 symmetry. The 
Hamiltonian in momentum space reads 

h= J2 *L^k*k, n k = ni®i2x2+n^°®a z , (5) 

kelBZ 

where * k = (*k,t>*k,4.) with ^ k;(T = 

(■SA,k,cr, Ssko-i sc.k.ff), an d where the 3x3 matrices 
and W]^ are given by 

/ -2tc x -2tCy \ 
H° k = -2tc x -tt'c x c y (6) 

\-2tCy -4t'C X Cy J 

and 

/0 \ 

•H k so = 4iAi SO -s x s v , (7) 

\0 S X Sy J 

with c p = cos(fc M /2) and s M = sin(fc M /2) (ji — x,y). The 
Hamiltonian matrix "H k consists of two uncoupled blocks cor- 
responding to the spin up and spin down projections, related 
by time-reversal symmetry, i.e., = (7iL k )* . Due to the 
time-reversal and inversion symmetries, the ISO coupling is 
unable to lift the spin degeneracy. For the calculation of the 
energy spectrum, it thus suffices to restrict our attention to 
one spin component, while keeping in mind that the resulting 
bands are doubly degenerate. In the following, we will restrict 
ourselves to the spin-up part of the Hamiltonian. We note that 
this Hamiltonian restricted to a single spin component is for- 
mally equivalent to the three-band model presented in Ref.[T5l 
for a copper oxide system. 

Let us first consider the case where only NN hopping is 
included (t' — Aiso = 0). The spectrum then consists of 
three bands (per spin component), where two dispersing bands 
touch the flat middle band via a cone-like dispersion at the M 
point, see Fig.|2ja). There is only one such point in the 1BZ. 16 
Since the Nielsen-Ninomiya theorem does not apply here due 
to the flat band, there is no fermion doubling. 17 When the ISO 
coupling is included (t' = and Aiso 7^ 0), a gap opens at 
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M, see Fig. gb). As shown in Refs. [18] and ]19\ the result- 
ing gaps are nontrivial, i.e., they exhibit helical edge states, 
and the upper and lower bands have nonzero Chern num- 
bers. These Chern numbers may be evaluated analytically, 
as demonstrated in Ref. Q3] The result of this computation 
is C T = (Cl,Cl,Cl) = (1,0,-1), for the Chern numbers 
of the lower, middle, and upper band, respectively, for the 
spin-up component. For the spin-down component, we have 
= —C'. From this property, we immediately obtain the 
Chern numbers for Aiso < from those for Aiso > by 
interchanging the roles for spin up and spin down. Since the 
Hamiltonian is diagonal in spin space, we can define the spin 
Chern number as C sp ; n = — C^. 20 Due to the symmetry 
between the two spin components, it holds that C sp ; n = 2C^. 
By virtue of the bulk-boundary correspondence, the spin Hall 
conductivity (in units of the spin conductivity quantum e/Air) 
equals 



0"SH 



= E ft. 

a: e Q <e F 



pin, a i 



(8) 



where the summation is over the filled bands (eF is the Fermi 
energy). The Z 2 -index v (defined as for an even and 1 for 
an odd number of edge-state pairs) is then related to spin Hall 
conductivity by v = ctsh/2 mod 2. An alternative way to 
derive the Chern numbers and the spin Hall conductivities is 
by diagonalising the system in a ribbon geometry, and to count 
the number of edge states that appear inside the bulk gapsP^El 
We remark that for systems that respect inversion symmetry, 
the Z2 invariant can also be calculated using the parity eigen- 
values at the time-reversal symmetric momenta. 22 The advan- 
tage of computing the spin Chern numbers as we do here, is 
that they are directly related to the quantized spin Hall con- 
ductivity in the band gaps. 

Now, let us investigate the effects of a purely real NNN 
hopping (£' 7^ and Aiso = 0). Unlike the ISO coupling, 
the real NNN term breaks the particle-hole symmetry of the 
spectrum. This case is, in some sense, more complicated than 
the previous one, because the Hamiltonian matrix cannot be 
written in the form R • L, with L a three-component vector of 
matrices that form a basis for a representation of su(2). How- 
ever, it is still possible to extend the components of L to a 
basis of su(3). This means that the Chern number does not 
have a simple interpretation as a winding number any longer. 
Moreover, the expressions for the eigenvalues are more com- 
plicated, although one can still find them analytically. 

We will now focus on the qualitative behavior of the dis- 
persions as function of t'. Let us start with t' = 0, for which 
the middle band is flat, as was already discussed before [see 
Fig.[3|a), dashed curves]. For t' > 0, the middle band starts 
to develop a maximum at the T-point [see Fig. [3|a), solid 
curves], and ultimately, this maximum touches the upper band 
for t'/t = 0.5 [Fig. [3}b)]. For t' /t > 0.5, the middle and 
upper bands touch at four inequivalent points in the 1BZ of 
the form Q = (±q,±q), where q is given by the condition 
cos(qa/2) = t/2t' [see Fig.[3jc)]. Similarly, the middle and 
lower bands touch for t' / 1 < —0.5, as shown in Fig.[3|d). The 
dispersion around these touching points Q [see FigT3|e)] re- 
sembles a tilted anisotropic Dirac cone, as shown in Fig.|3jf). 
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FIG. 3. (Color online) Dispersions for various values of t' along the 
high symmetry lines in the 1BZ. (a) The case where t' = is shown 
by the dashed curves. For < i! < 0.5t the middle band starts 
to develop a maximum at the T-point (solid curves), (b) For t' = 
0.5 the maximum at the F-point touches the upper band, (c) When 
t' > 0.5t a tilted anisotropic Dirac cone forms at the Q-point. (d) 
The dispersion for t' < — 0.5t. In (e), the four corners of the dotted 
square indicate the positions of the four tilted anisotropic Dirac cones 
in the 1BZ. (f) A three-dimensional picture of the anisotropic tilted 
Dirac cone at Q. 



The four solutions reflect the C4 symmetry of the underlying 
lattice and the fact that they do not coincide with high symme- 
try points in the 1BZ. Moreover, we observe that the position 
of these cones in the 1BZ is tunable by t'. In the vicinity of 
the touching points, the dispersion can be approximated by 
the generalized Weyl HamiltoniarPl 



#eff(p) = E t 2x2 +t(^p 



't 



1 



1 



2x2 



-OyV 



y/3 ' 3' 

(9) 

where p ± — p x ±p y withp = (p x ,p y ) = k—Q is the momen- 
tum with respect to the touching point Q,t = t-i/l — (t/2t') 2 , 
and Eq = t/t' is the energy of the bands at Q. This expression 
may be found using the Luttinger-Kohn representation 2 ^of the 
bands at the band-touching points. The tilt comes from the 
term proportional to p + l 2X 2- We remark that the maximum- 
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FIG. 4. (Color online) (a) The black dashed curves indicate the dis- 
persions for Aiso = and the blue continuous curves indicate the 
situation for Aiso 7^ 0. There is no full gap when Aiso is below 
a critical value A[ S0 . (b) For a sufficiently large Aiso > A[ SO , the 
spectrum is fully gapped. 



tilt condition as defined in Ref. [23] is not satisfied for these 
cones. In addition, the Hamiltonian Q is only of limited in- 
terest for our purposes because by itself it does not reproduce 
the topological properties that we analyze. 

The most general case is when both the real and imaginary 
parts of the NNN hopping are included, t! ^ and Aiso 0- 
In this case, the ISO term opens a gap at M. Regarding the 
Chern numbers of the bands, we immediately understand the 
case t' < 0.5t: Such a state is adiabatically connected to the 
case for which t' = 0, and no gaps close in this regime, so 
that the Chern numbers of the bands are exactly the same as 
if one would put t' = 0. The situation is however more subtle 
for t! > 0.5t. As mentioned previously, a band touching point 
occurs at the T point for Aiso = 0. Interestingly, this still oc- 
curs in the presence of ISO coupling, since its contribution to 
the Hamiltonian [see Eq. Q] vanishes at V. Thus, the T point 
is protected from the effects of the ISO coupling. Because the 
upper and middle band do touch when t! is increased to 0.5t, 
the Chern numbers of the bands may change. An analytical 
evaluation of the Chern numbers is in principle possible, but 
the expressions become tedious to deal with. For this reason, 
we choose to compute the Chern numbers numerically by re- 
sorting to the method proposed by Fukui et a/P^For the case 
t! = 0, the analytical results from a computation similar to the 
one presented in Ref. 15 coincide with our numerical findings. 

Before we present the results from this calculation, we 
would like to point out a subtlety that arises in the interpre- 
tation of the Chern numbers in terms of the spin Hall con- 
ductivity. For the latter quantity to be quantized, there must 
be a full gap in the spectrum, i.e., there must be a range of 
Fermi energies for which there are only edge states and no 
bulk states. However, as shown in Fig. |4f a), there may be no 
full band gap between the middle and upper bands, although 
they are separated at each point in the 1BZ. In this situation, 
i.e., if the minimum of the upper band lies at a lower energy 
value than the maximum of the middle band at a different mo- 
mentum (eo(r) > e + (M)), we say that the spectrum exhibits 
a negative indirect gap. The Chern numbers of the bands are 
well-defined, since they do not touch anywhere, and no topo- 
logical phase transition occurs. In this example, a 2/3-filled 
system will always have both the upper band and the mid- 
dle band partially filled, which would classify the bulk as a 
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FIG. 5. (Color online) Phase diagram as function of Aiso and t' . 
The red (solid) lines indicate positions for which a band gap closes, 
meaning that the Chern numbers can change. In the regions enclosed 
by these lines the various (spin) Chern numbers = C S pm/2 are 
indicated. The regions delimited by the blue (light gray) or green 
(dark gray) dashed lines indicate the parameter regimes for which 
there is no full gap between the middle and upper bands or middle 
and lower bands, respectively. 



semimetal, preventing the helical edge states to be observed. 
In the semimetallic regime, i.e., with partially filled bands, the 
spin Hall conductivity is not quantized.^ For 1/3-filling, this 
problem does not occur since in this case there is a full band 
gap. At this filling, the system behaves as an insulator and the 
spin Hall conductivity carried by the helical edge states can 
be experimentally observed. By virtue of symmetry, a similar 
reasoning occurs for t' < 0, but with the roles of the 1 /3 and 
2/3 filling, and of the lower and upper bands interchanged. 
Finally, the spectrum is fully gapped if the ISO coupling is 
sufficiently strong, i.e., there is a critical value Aj SO (i'), such 
that the system exhibits only full gaps if the ISO coupling sat- 
isfies I Aiso I > Af so , see 

Fig.gb). 

In Fig. [5] we display the phase diagram in the (t' , Aiso)- 
plane, where we indicate the regions with distinct Chern num- 
bers for the bands. The topological phase transitions, i.e., 
the values (t' , Aiso) where band touchings occur, are indi- 
cated by solid lines. Moreover, we have shaded the parameter 
regimes in blue (green) where the middle and upper (lower) 
bands overlap in energy and the only full gap is at 1 /3 (2 /3) 
filling. In all open gaps, the system behaves as a quantum spin 
Hall insulator with spin Hall conductivity equal to osh = ±1- 
For \t' /t\ < 0.5, the conductivities in the lower (1/3 filling) 
and upper (2/3 filling) gap are equal, since the middle band 
has zero Chern number. For \t' jt\ > 0.5, the two bands have 
opposite conductivities. 

Every phase transition in the phase diagram may be under- 
stood in terms of the difference in the Chern numbers of the 
phases at both sides of the transition. For the transitions de- 
fined by the lines t'/t = ±0.5, either the lower or the up- 
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per gap closes, and the system behaves as a metal at the gap 
closing energy; the other gap is still helical. For example, if 
if /t = 0.5 and Aiso > 0, the difference of the Chern numbers 
is AC r = (1, 0, -1) - (1, -2, 1) = (0, 2, -2). We observe 
that the change of the Chern numbers of the two touching 
bands is ±2. This value can be understood from the fact that 
the bands touch at T, and they behave quadratically around 
this point, i.e., there is a Berry phase of 4tt associated to this 
touching point. For the transitions at Aiso = (i.e., between 
the topological phases for Aiso > an d Aiso < 0), we dis- 
tinguish two cases. First, for \t' /t\ < 0.5, the three bands 
touch at M, and the Chern-number difference between the two 
phases is equal to AC T = (1, 0, -l)-(-l, 0, 1) = (2,0,-2), 
i.e., ±2 for the upper and lower bands, respectively. Secondly, 
for t' /t > 0.5, the three bands touch at M like before, but also 
at the four Weyl cones. The change in Chern numbers is there- 
fore AC T = (1, -2, 1) - (-1, 2, -1) = (2, -4, 2), with the 
contributions (2, 0, —2) from the M point and (0, —4, 4) from 
the four Weyl points, respectively. The transition at Aiso = 
for if jt < —0.5 is similar. 



B. The Lieb lattice with a dimerization term 

In Ref. [T8l it has been shown that a dimerization term, de- 
fined by changing the hopping alternatingly to t ± a, leads to 
a trivial gap. In this section, we investigate the effects of the 
dimerization term in combination with the other terms in the 
Hamiltonian ([l}. The Hamiltonian then reads 



u^ = ul® i 2x2 + u l £° ®a z + ni im 

where the dimerization term is given by 



tlx. 




—2ias T —lias 



(10) 



(11) 



in momentum space. This term breaks the C4 symmetry of 
the lattice. In principle, the signs before s y may be chosen 
oppositely, but this choice does not affect the following results 
qualitatively, since it amounts merely to the transformation 
(k x , k y ) — > (k x , —k y ) in momentum space. 

The effect of the dimerization for a/t = 0.3 on the phase 
diagram is shown in Fig. [61a). We observe that, compared to 
the case for a = in Fig. [51 there are additional phase tran- 
sitions; in particular, the existing phase transitions have been 
split up. Comparing Figs. [5] and [6| a), we observe that in both 
cases there are phase transitions at t' ji = ±0.5. However, for 
a^0, there is an additional phase transition close to the one 
at if /t = ±0.5 (for |Aiso/*| > 0.2). Let us compare the dif- 
ference in Chern numbers at these transitions as before: For 
a = 0, this difference is AC T = (1,0,-1) - (1,-2,1) = 
(0, 2, —2). With nonzero a, there are two transitions, from 
(1, 0, -1) to (1, -1, 0) and from (1, -1, 0) to (1, -2, 1); thus, 
the dimerization term has created an intermediate phase with 
= (1,-1,0), which has a helical structure in the lower gap 
(unaffected by the transition), and a trivial upper gap. The 
first transition [from (1,0,-1) to (1,-1,0) at \t' /t\ = 0.5] 
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FIG. 6. (Color online) (a) Phase diagram as function of Aiso and t' 
for a nonzero dimerization term (a/t = 0.3). We indicate the phase 
transitions by red (solid) lines, and the Chern numbers of the bands 
— Cspin/2 are indicated in the regions delimited by these lines. In 
the shaded regions, one gap is helical and the other is trivial, (b) The 
regimes with indirect negative gaps (semimetallic behavior) at 1/3 
or 2/3 filling are indicated by the green (dark gray) and blue (light 
gray) regions. In the white regions, the system has full band gaps at 
1/3 and 2/3 filling. 



occurs due to a band touching at the point T; expanding the 
bands around this point, we observe that they are quadratic 
along one diagonal [k = (q, q)] and linear along the other 
[k = (q, — 9)]. The other phase transition [from (1,-1,0) 
to (1, —2, 1) for t' )t slightly larger than 0.5] occurs due to a 
band touching involving a single Weyl cone on the diagonal 
k = (q, q). Thus, the original transition with a Chern-number 
difference of ±2 for the touching bands has been split up in 
two transitions with a Chern-number difference of ±1. 

The transition at zero ISO coupling, originally between 
the phase with (1,0,-1) and (—1,0,1) has also been split 
up by the dimerization term. Going from Aiso = 0.3 to 
A ISO = -0.3 for < \if/t\ < 0.5 and a/t = 0.3 , one 
encounters four gap closures (and two for if = 0), where the 
bands (the middle and either the upper or lower) touch in one 
point (in a Weyl cone-like behavior). At Aiso = there is 
no transition, i.e., the system remains trivially gapped. For 
t! /t > 0.5, the dimerization splits the transition originally at 
Aiso = into three: For a > 0, the middle and upper bands 
touch at two Weyl cone-like points if Aiso = 0, and the corre- 
sponding change of Chern numbers is ±2. At Aiso A ~ ±0.2 
(for a/t — 0.3), the middle and lower bands touch at one 
point, accounting for a change of Chern numbers of ±1. 

The dimerization term also modifies the parameter regimes 
where one of the gaps is a negative indirect gap, as may be 
observed by comparing the shaded regions of Figs. [5] and 
Fig. |6|b) (without and with dimerization term, respectively). 
There are two qualitative differences: First, with the dimer- 
ization term, the critical spin-orbit coupling Aj go (i')' tends 
to infinity for \t'/t\ —> 0.5, whereas Aj SO (t') remains fi- 
nite if a = 0. Secondly, with dimerization, there are no 
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negative indirect gaps if the NNN hopping amplitude satis- 
fies \t'\ < H/%/2 (i.e., \t'/t\ < 0.21, approximately, for 
a/t = 0.3). We remark that for all topological phases in- 
dicated in Fig. |6ja), there is a regime of parameters (1f, Aiso) 
for which it may be observed in a system with two full band 
gaps, i.e., such that the spin Hall conductivity is quantized in 
both gaps. 



III. THE KAGOME LATTICE 

The existence of topological phase transitions driven by 
the NNN hopping t' in the Lieb lattice raises the question 
whether they are specific to that lattice or whether they would 
also appear in other lattices. In this section, we investigate 
the kagome lattice, which has the lattice structure shown in 
Fig. |7ja). The unit cell consists of three sites, but these all 
have the same number of NNs, unlike the Lieb lattice. In mo- 
mentum space [for the structure of the 1BZ, see Fig.|7jb)], the 
Hamiltonian for this system is given by 



Hi = <- NN ®l 



in the basis 



*7 /K,NNN ~ -, 



-H 



K,ISO„ 



(12) 



(*k,t,*k,4.). where * k ,, 



f - - \ o o 1 /K,NN n ,K,NNN 

(s J 4,k,ff,Sfl,k.j,sc.k,ff)- The 3 x 3 matrices % k ' ,« k ' 
and rt k are given by 



t,K,NN _ 



n 



K,NNN 



(0 ci c 3 N 
-2t ci c 2 
\C3 c 2 

f0 bi b 3 ■ 
-2t' [h b 2 
\h b 2 



(13) 



(14) 




K M 



FIG. 7. (Color online) (a) The kagome lattice with the three inequiv- 
alent sites A, B, and C inside the unit cell (dashed area). The NN 
and NNN hoppings are indicated by the red arrows, (b) First Bril- 
louin zone of the kagome lattice with high-symmetry points, (c) Dis- 
persions along contour TKMY for Aiso = and t' = (dashed 
curves) and for Aiso A = 0.1 and t' = (solid curves). (d,e,f) 
Dispersions for Aiso/t = 0.1 and t'/t = 0.39, t'/t = 0.443, and 
t'/t — 0.6, respectively. The upper gap closes for t'/t ~ 0.443. 



and 






bi 




2«Aiso ~h 





h 


V ^3 


-b 2 






where Cj = cos(k • e,) and bi = cos(k • f,) are given in 
terms of the NN vectors ei = (1,0), e 2 = (-1, \/3)/2 and 
e 3 = (— 1, — \/3)/2, and the NNN vectors fi = e 2 - e 3 , 
f2 = e 3 — ei, and f 3 = ei — e 2 . 

In absence of NNN hopping and ISO coupling, the spec- 
trum consists of three bands: Two dispersing bands, which 
resemble the dispersion of graphene, with two inequivalent 
Dirac points inside the 1BZ, where the bands touch at E/t = 
— 1, and a flat band at E/t = 2, which touches the middle 
band at rj23see Fig.|7jc), dashed lines. If only real NNN hop- 
ping is included (i.e., t' ^ and Aiso = 0), then no gaps open 
between the bands. If only ISO coupling is included [see the 
solid lines in Fig. |7jc)], then gaps exhibiting helical states ap- 
pear. Unlike the Lieb lattice, varying the ISO coupling drives 
several phase transitions, not just one for Aiso = 0- For the 
case that both the real and imaginary parts of the NNN hop- 
ping are nonzero (t' ^ and Aiso 7^ 0), we find many band 



touching points, see Figs. |7Jd,e,f). For instance, one of the 
gaps closes at K and K' for Aiso/* = ±V3(*7* - 1/2). Ad- 
ditional gap transitions occur due to bands touching at points 
between T and K, between T and M, and between M and 
R( >. Without explicit computation of the Chern numbers, 
we can use the symmetry properties of the 1BZ to predict the 
change of Chern number in these transitions. The band touch- 
ing points at T, at \ between T and K, between T and 
M, and between M and ' have multiplicities 1, 2, 6, 6, 
and 6, respectively. The difference in the Chern numbers of 
the bands between the phases at both sides of the correspond- 
ing transitions is an integer multiple of this multiplicity. Al- 
though a phase diagram like Fig. [5] can be made also for the 
kagome lattice, it is outside the scope of this paper. The ques- 
tion whether topological phase transitions driven by real NNN 
appear in other two-dimensional Dirac-like systems than the 
Lieb lattice has been answered positively. Next, we investi- 
gate the honeycomb lattice, aiming at checking how generic 
this behavior is. 
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IV. THE HONEYCOMB LATTICE 

The honeycomb lattice [see Figs. |8ja) and (b)] is one of 
the most intensively studied Dirac systems because it models 
the electronic structure of graphene very well. Since the work 
by Kane and Melep^l it j s known that ISO coupling in the 
honeycomb lattice gives rise to the quantum spin Hall phase. 
The real NNN hopping has been neglected in this model, al- 
though it is several orders of magnitude stronger than the ISO 
coupling: In pure graphene, \t' /t\ ~ 0.1 and \\so/t ~ 10 -6 - 
10 _4 P^Here, our aim is to study the dispersions and the pos- 
sibilities for phase transitions is both ISO coupling and real 
NNN hopping are included in the tight-binding model for the 
honeycomb lattice. 

The Hamiltonian in momentum space for the honeycomb 
lattice is given by 



O/H _ T/H.NN 



12x2" 



-H 



H.NNN , 



12x2 



+H 



HJSO, 



in the basis = (^k,-^ ^k,j.), where 

(sA.k,cr, SBM.a)- The 2x2 matrices 7^' NN , H k 
7i k are given by 



u 



u 



and 



H.NN 
k 

H,NNN 



s(k) 
5*(k) 

h(k) + h*{k) 






h(k) + h*(k) 



(16) 
, and 

(17) 
(18) 



H 



HJSO 



= -iX 



ISO 



fh(k) - h*(k) 





-h(k) + h*(k) 



(19) 



where we define g(k) = e 



-k-ei 



Q ike 2 



o*ke 3 



oikfi 



oikf 2 



oik- f a 



and h(k) 



, with ej and the same as given for the 
kagome lattice. For this model, it is straightforward to derive 
the dispersion of the two bands analytically, as 



E = t'[h(k) + h*(k)]±^\tg(k)\* + X 2 ISO [ih(k) - ih*(k)}2. 

(20) 

The bands touch if and only if the square root equals zero. 
Since both terms in this square root are positive, the bands 
touch only when g(k) = (i.e., at k = K or k = K') and 
AisoIm/i(k) = 0; because lmh(K) ^ 7^ Imh(K'), the 
bands can only touch for Aiso = 0. The real NNN hopping 
has no effect on the touching of the bands; it merely shifts the 
energies of both bands up or down. Thus, the conclusion is 
that in the honeycomb lattice, there are no phase transitions 
driven by the real NNN coupling. The only topological phase 
transition in this model is the one at Aiso = 0, which forms 
the boundary between two helical phases with opposite signs 
of the spin Hall conductivity. 

However, for the honeycomb lattice in a magnetic field, the 
results are different. (For the framework of a honeycomb lat- 
tice in a magnetic field, see e.g. Refs. l9land[T0l) The Chern 
numbers and are well-defined because of the absence 
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E/t 
1.9 
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S ~ 













FIG. 8. (Color online) (a) The honeycomb lattice with the two in- 
equivalent sites A and B inside the unit cell (dashed area). The 
NN and NNN hoppings are indicated by the red arrows, (b) First 
Brillouin zone of the honeycomb lattice with high-symmetry points, 
(c) Dispersions for Aiso = 0, t' = 0, and 1/3 flux quantum per 
unit cell, along the contour TKMT of the magnetic Brillouin zone 
(1/3 of the size of the 1BZ). (d,e,f) Edge-state plots in the energy 
regime of the two highest bands for Aiso /£ = and t' /t = —0.18, 
t'/t = —0.2, and t' jt = —0.22, respectively. For these plots, the 
dispersions have been computed for the system in a cylindrical ge- 
ometry (see Ref.llOt. The light gray curves are the bulk bands. The 
red (dark gray) and blue (medium gray) curves are the edge states on 
the two opposite edges of the cylinder. The variation of t' closes the 
gap between these bands, and changes the Hall conductivity in this 
gap from — 2e 2 /h to 4e 2 fh. 



of a coupling between the spin components, but they do not 
satisfy the condition = — C^, which applies only if there is 
time-reversal symmetry. 

In the case that ISO coupling is absent, there are no terms 
that break the spin symmetry, so that one observes spin- 
degenerate chiral (i.e., quantum Hall) phases only. Due to 
the spin degeneracy, the Chern numbers satisfy = C , 
and the Hall conductivity is quantized in steps of 2e 2 jh. In 
Fig. |8jc), we have displayed the band structure for the hon- 
eycomb lattice under a flux of 1 /3 flux quantum per unit cell. 
In this case, there are 6 bands per spin component, rather than 
2. Here, in absence of ISO coupling, the variation of the real 
NNN hopping induces phase transitions between chiral states. 
In Figs.[8|d,e,f), we show that the gap between the two high- 
est bands closes when varying t'/t. In Figs.[8jd) and (f), there 
are 1 and 2 edge-state pairs in the relevant gap, respectively. 
In addition, the direction of the edge states is flipped. Thus, 
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the Hall conductivity inside this gap changes from uh = —2 
to oh =4 (in units of e 2 /h; recall that the spin degener- 
acy yields an additional factor of 2). Here, we note that the 
change of 6 is explained from the fact that there are 3 inequiv- 
alent band touching points in the 1BZ due to the fact that there 
is 1/3 flux quantum per unit cell, and from the two-fold spin 
degeneracy. 

In presence of ISO coupling together with the magnetic 
field, more exotic phase transitions can be expected. In that 
case, the spin degeneracy is lifted and the condition = 
±C^ is generally not satisfied. As a consequence, more types 
of topological phases can be observed than the helical and 
spin-degenerate chiral phases, e.g., chiral phases which are 
non-spin degenerate, and even phases which are neither chiral 
nor helical, but a mixture of both. In Refs. [9] and [TO] it has 
been shown that variation of the ISO coupling leads to a topo- 
logical phase transition between a chiral and a helical state, 
as well as several transitions between other types of topolog- 
ical phases. If additionally real NNN hopping is present, this 
picture remains qualitatively unaltered. In conclusion, topo- 
logical phase transitions in the honeycomb lattice subjected to 
a perpendicular magnetic field can be driven by ISO coupling 
as well as by real NNN hopping. 

V. DISCUSSION 

The theoretical results presented above raise the question 
whether the transitions driven by real NNN hopping appear 
in an experimentally accessible regime. We envisage sys- 
tems with cold atoms in an optical lattice to be the most 
promising candidates for observation of these phase transi- 
tions, because of the possibility to tune the hopping parame- 
ters. In a static optical lattice, t' /t is typically of the order of 
10 _2 P^far from the transitions that happen at t'/t = 0.5. 
However, shaken optical lattices exhibit renormalized hop- 
ping coefficients, where the bare values become multiplied 
by a Bessel function. 1 ^ In the most studied case, one applies 
a periodic shaking to the lattice with the effective potential 
V = Vo cos(u;t) Yli,j(^ + j)fkj> where Vq is the ampli- 
tude and u is the frequency of the shaking. Here, la- 
bels the lattice sites in the two-dimensional lattice, r is real 
time, and denotes the density operator at site In 
the high-frequency limit hui Vq, t, the effective hopping 
may be determined by applying Floquet theory. One finds that 
t c s = tJo(K), where the shaking parameter K = Vo/tiw. 
Moreover, the argument of the Bessel function generally dif- 
fers for NN and NNN hopping (see e.g. Ref.[l2|for the case of 
a non-separable square lattice), thus allowing one to tune the 
ratio t' cS /t c g in the entire range of possibilities (0 to oo). 

The detection of topological phases in optical-lattice exper- 
iments presents a challenge because it is impossible to per- 



form (Hall) conductivity measurements in such systems. First 
of all, the atoms are neutral, so that there are no electric cur- 
rents. Therefore, different detection methods must be envis- 
aged to image the edge states and to determine the topolog- 
ical properties. Recently, there has been a proposal to de- 
tect chiral edge states in a Hofstadter optical lattice by using 
angular- momentum sensitive Bragg spectroscopy. 30 Addition- 
ally, the edge states are probed in a clever way using Raman 
transitions, in order to isolate the edge-state signal from the 
background. A second potential problem is the harmonic trap 
in which the atoms are confined. The lack of sharp edges 
means that the notion of edge states is no longer well-defined. 
This phenomenon is visible in the spectrum of a system in a 
harmonic background potential: The ordered structure of the 
spectrum, consisting of the bulk bands and a few edge states in 
the gaps, is destroyed by this potential. 31 However, it has been 
shown that the topological invariants remain well-defined in 
the case of a system with soft edges. In other words, sharp 
edges are not necessary in order to observe the topological 
phases.™ 

Besides the more complex optical lattices with multi-site 
unit cells considered here, interesting topological phases may 
also occur in simpler square lattices loaded with bosons in 
higher excited states, e.g., p, d, and / orbitals (see Refs. [32] 
1331 and [34] respectively). Indeed, for a square superlattice 
with deeper and shallower wells arranged on a checkerboard 
pattern, it is possible to realize relatively long-lived Bose- 
Einstein condensates in the deeper wells which are locally in 
a metastable p x or p y state, whereas the ones in the shallower 
wells are in an s-orbital state. The single-particle spectrum 
of this system is then very similar to the one of the Lieb lat- 
tice, which is understood from the formal equivalence of the 
Hamiltonians involving three orbitals per unit cell. 32 

In conclusion, we have shown that in several Dirac-type 
tight-binding models, it is possible to drive topological phase 
transitions with the real NNN hopping. In the Lieb and 
kagome lattices, it drives transitions between helical phases 
with different values of the spin Hall conductivity. This is 
not true for the honeycomb lattice in zero magnetic field. In 
a magnetic field, the phase transitions are possible. In this 
aspect, we conjecture that it is necessary to have at least 
two gaps (three bands) in order for these phase transitions to 
appear. Our results could be observed with state-of-the-art 
shaken optical lattices loaded with fermions. 
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